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Abstract 

In this work we have used extensive Monte Carlo calculations to study the planar to paramag- 
netic phase transition in the two-dimensional anisotropic Heisenberg model with dipolar interac- 
tions (AHd) considering the true long-range character of the dipolar interactions by means of the 
Ewald summation. Our results are consistent with an order-disorder phase transition with unusual 
critical exponents in agreement with our previous results for the Planar Rotator model with dipolar 
interactions. Nevertheless, our results disagrees with the Renormalization Group results of Maier 
and Schwabl [PRB, 70, 134430 (2004)] and the results of Rapini et. al. [PRB, 75, 014425 (2007)], 
where the AHd was studied using a cut-off in the evaluation of the dipolar interactions. We argue 
that besides the long-range character of dipolar interactions their anisotropic character may have 
a deeper effect in the system than previously believed. Besides, our results shows that the use of a 
cut-off radius in the evaluation of dipolar interactions must be avoided when analyzing the critical 
behavior of magnetic systems, since it may lead to erroneous results. 

PACS numbers: 75.40.Cx, 75.40.Mg, 75.10.Hk 
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I. INTRODUCTION 



Magnetism is one of the most studied subjects in physics. In particular, the study of 
the statistical behavior of spin models was the trend in the years 1960's and 1970's. They 
gave fruitful contribution to the understanding of several phenomena, not only in Physics, 
but in several other fields. The concept of symmetry breaking and universality transcend 
the subject of physics to many other areas, being thus of great importance. More recently, 
mainly due to the growing interest in magnetic thin-films, magnetic nano-dots and arrays 
of magnetic nanoparticles -~— , there were a renewed interest in the study of spin models. 

In many magnetic systems a long range dipole-dipole energy term has to be considered 
beside the exchange interaction between neighboring sites and the anisotropics present in the 
system. The study of such models is mainly associated with the development of magnetic- 
nonmagnetic multilayer for the purpose of giant magnetoresistence applications and mag- 
netic nanoparticle arrays as an alternative for new media storage. In addition, experiments 
on epitaxial magnetic layers have shown that a huge variety of complex structures can de- 
velop in the system^^. Rich magnetic domain structures like stripes, chevrons, labyrinths, 
and bubbles associated with the competition between dipolar long-range interactions and 
a strong anisotropy perpendicular to the plane of the film were observed experimentally. 
A lot of theoretical work has been done on the morphology and stability of these magnetic 
structures^— . Beside that, it has been observed the existence of a switching transition from 
perpendicular to in-plane ordering at low but finite temperature^: at low temperature the 
film magnetization is perpendicular to the film surface; rising temperature the magnetization 
flips to an in-plane configuration. Eventually the out-of-plane and the in-plane magnetiza- 
tion become zero. Although the structures developed in the system are well known, the 
phase diagram of the model is still not completely understood. 

Here, our interest is in magnetic thin-films with an out-of-plane anisotropy as described 
above. For such a system we can write a model Hamiltonian as follows 



H = - J E $ • s s) - A E(^) 2 + D E 
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known as anisotropic Heisenberg model with dipolar interactions (AHd). Here, we consider 
a ferromagnetic system, so that J > 0, A is an easy-axis anisotropy, and D is the strength of 



2 



1 

0.8 
0.6 

T 

0.4 
0.2 

0.1 0.12 0.14 0.16 0.18 0.2 
D/J 

Figure 1: (Color online) Phase diagram of the anisotropic Heisenberg model with dipolar interac- 
tions (AHd) for fixed A/J = 1 in the (D/J,T) space. The black solid line represents the transition 
lines as obtained using a cut-off in the dipolar interactions^ and the solid black line are the results 
obtained when full long-range interactions are considered by means of the Ewald summation (this 
work). The phase I is an Ising-like phase characterized by an ordered out-of-plane alignment of 
spins (that may present stripe- like configurations for full long-rang interactions). Phase II is an 
ordered planar ferromagnetic state and phase III is a paramagnetic one. 




the dipolar interactions, r^- is a vector connecting sites i and j while < i, j > means that the 
first summation is to be evaluated for nearest neighbors only. For the dipolar interactions 
the summation is evaluated over all pairs i j and in the single-ion anisotropy term the 
summation is evaluated for all sites in the lattice. For D ^ the system is frustrated 
due to the competition between the dipolar and the anisotropic terms. For small D/J 
compared to A/J we can expect the system to have an Ising-like behavior since an out-of- 
plane configuration of the spins is expected. If D is not too small we can expect a transition 
of the spins from out-of-plane to in-plane configuration. For large enough D out-of-plane 
configurations become unstable such that, the system lowers its energy by turning the spins 
into an in-plane ferromagnetic arrangement. Earlier works on this model, which discuss 
the phase diagram, were mostly done using renormalization group approach and numerical 
Monte Carlo simulation^ - — . They agree between themselves in the main features. The 
phase diagram for fixed A/J is schematically shown in Figfl]in the (D/J,T) space. From 
Monte Carlo (MC) results it is found that there are three regions labeled in Fig. [1] as I, 
II, and III. Phase I corresponds to an out-of-plane magnetization, phase II has in-plane 
magnetization, and phase III is paramagnetic. The border line between phase I and phase II 
is believed to be of first order and that between regions I and II is a second order one^. The 
transition line between regions II and III has a non clear character. Some authors reported 
they found a second order line^. In reference 
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the author's claim that the transition is of 
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the BKT type. In many cases a second order phase transition can be confused with a BKT 
transition. This can be due to the analysis of the finite size effects to be done in not large 
enough lattices. Although the different results point out in the direction of a second order or 
a BKT transitions^- between region II and III, much care has to be taken because they were 
obtained by using a cutoff radius r c in the dipolar interaction. The long-range character of 
the potential is lost and then long-range order may not develop^. As a consequence, it will 
not be surprising if a completely different scenario emerges when considering the true long- 
range dipolar interaction, since long-range order is expected to be present^. In a recent 
paper, Maier and Schwab!^ analyzed the phase transition in the dipolar planar rotator 
(dPR) model via renormalization group techniques. It is expected that the dPR model 
describes the critical behavior of the planar to paramagnetic transition in the anisotropic 
Heisenberg model with dipolar interactions (AHd) since they have the same simetry. Their 
results indicate that the dPR model belongs to a new universality class characterized by 
an exponential behavior of the magnetization, susceptibility and correlation length. Besides 
that, the specific heat was found to be non- divergent, as occurs in the BKT phase transition. 
More recently M61 and Costa^ 1 ^ studied, by using Monte Carlo simulations, the AHd model 
in a bilayer system and the dPR model. In both cases they found strong evidences for the 
transition between region II and III to be in another universality class characterized by 
a mixed behavior between order disorder and BKT transitions. Besides, even in a 2D 
dipole lattice the scenario is not clear due to some conflicting results (see, for example, 
references^— ) . It should be noticed however that in most of these works the full long-range 
character of dipolar interactions are not taken into account properly. Indeed, the dipolar 
interaction is conditionally convergent in two dimensions due to its anisotropic character, 
such that the use of the minimal-image convention or the introduction of cut-offs in the 
potential may not describe the real behavior of the system. 

In this work we have used extensive Monte Carlo simulations to study the transition line 
between regions II and III (the planar to paramagnetic phase transition) in the easy-axis 
Anisotropic Heisenberg model with dipolar interactions (AHd) (see EqnJTJ). Although the 
system in consideration has a significant out-of-plane anisotropy our main interest is in a 
easy-plane symmetry brought about by the dipolar interactions, so that for a suitable choice 
of the model parameters the spins lie in the film plane. Indeed, our goal is to compare 
the results obtained for the above mentioned model when full long-range dipolar interac- 
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tions are considered by using the Ewald summational 3 ^ (this work) with those where a 
cut-off radius is introduced in the evaluation of dipolar interactions^. Our results clearly 
indicate that considering the long-range character of dipolar interactions the transition line 
between the planar and paramagnetic phases is of the order-disorder type, being charac- 



exponents, instead of being a 
This observation can be of 



21 



terized by a non-divergent specific heat and unusual critica 
Berezinskii-Kosterlitz-Thouless one as reported in reference 
great importance since the effects of dipolar interactions is enhanced in systems in the sub- 
micrometer scale, as magnetic nano-dots^ 3 -^ 13 -^ 3 ^ and arrays of magnetic nanoparticles^ 3 ^ 
and can shed extra light on some controversial results in the literature^ - — . 



II. SIMULATION BACKGROUND 



26 



27. The reader is referred 



In this work we followed the same methodology used in Refs. 
to these works for a more detailed description, specially about the determination of critical 
temperature and exponents. Our Monte Carlo procedure consists of a simple Metropolis 
algorithm^! where one Monte Carlo step (MCS) consists of an attempt to assign a new 
random direction to each spin in the lattice. To equilibrate the system we have used 100 x L 2 
MCS which has been found to be sufficient to reach equilibrium, even in the vicinity of the 
transition. In our scheme, two sets of simulations have been performed. In the first one, 
we preliminarily explored the thermodynamic behavior of the model in order to estimate 
the position of the maxima of the specific heat and susceptibilities and the crossings of the 
fourth order Binder's cumulant. In this first approach we used lattice sizes in the interval 
20 < L < 50. Once the possible transition temperature is determined, we refined the results 
by using single and multiple histogram methods. We produced the histograms for each 
lattice size in the interval 20 < L < 120 and they were built at/close to the estimated 
critical temperatures corresponding to the maxima and/or crossing points obtained in step 
1. To construct the histograms at least 2 x 10 7 configurations were obtained using at least 3 
distinct runs. These histograms are summed so that we obtain a new histogram that allow 
us to explore a wider range of temperature (an example of the use of histograms can be 



found in Ref. l26l ). Periodic boundary conditions are assumed in the directions x and y. 
To take into account the long range character of the dipolar interaction we use the Ewald 
summation to calculate the energy of the system 3 ^ 3 - 3 -. 
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All simulations where done using a square lattice, A/ J = 1 and D = 0.3 J. Energy 
was measured in units of J and temperature in units of J/ks, where ks is the Boltzmann 
constant. Our choice of D = 0.3J was to guarantee that the planar behavior of the system 
was not much affected by the frustration existent near the multicritical point where the three 
lines shown in Fig. [T] come together. We have devoted our efforts to determine a number 
of thermodynamic quantities, namely the specific heat, magnetization, susceptibility, fourth 
order Binder's cumulant and moments of magnetization as described elsewhere^ 1 ^. 



III. SIMULATION RESULTS 



Concerning the systems' magnetization no significant size dependence is observed in low 
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where a cut-off radius were used in the 



temperatures, unlike the results shown in Ref. 
evaluation of dipolar interactions. This may be an evidence that as the full long-range char- 
acter of dipolar interactions are taken into account long-range order develops, as expected 
by the results of Maleev^. 

In figure [2] we show a log — log plot of the maxima of the susceptibility as a function of 
the lattice size for L = 20,40,80 and 120. The data are very well adjusted by a straight 
line with slope j/u = 1.763(1) exhibiting a power law behavior. This value of the exponent 
7/V is quite near the expected one for a transition in the Ising universality class (1.75). 
Considering the Ising universality class we were able to determine the critical temperature 
by using the location of the maxima of the specific heat and susceptibility and the crossing 
point of the Binder's cumulant, which gives T^ sing = 0.946(1). By using this value and 
plotting \tl(Mxy x hi(L) at T = T c we have found (3/v == 0.163(6), which is quite different 
from the expected value for the Ising universality class (0.125). 

The last result may indicate that the assumption of the Ising universality class may 
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, we 



not be correct. Indeed, by analyzing the moments of magnetization defined in Ref. 
obtain 1/u — 0.82(2) and T c j = 0.943(1). This value of the exponent u contrasts with the 
expected for the Ising universality class, although the value for the critical temperature is 
approximately the same. Reanalyzing our previous estimates for the critical temperature 
obtained using the location of the specific heat and maxima of susceptibilities using this 
new value of the exponent v we obtain: T c c " = 0.945(1) and T c x = 0.943(1) (it is worthy to 
note that in the analysis of the specific heat data the point corresponding to L = 20 was 
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Figure 2: (Color online) Log-log plot of the maxima of the planar susceptibility as a function 
of the lattice size. The solid red line shows the best linear fit of the data given the exponent 
7/1/ = 1.763(1). The error bars are shown inside the symbols. 



disregarded in both cases). Looking to the crossing point of the Binder's cumulant we have 
found T^F 4 = 0.944(2). We have thus, as our new estimate for the mean critical temperature 
T c = 0.944(1). Using this new value of the critical temperature we obtain /3/u = 0.149(7) 
in the analysis of the magnetization data. 

To distinguish between these scenarios in figure |3] we show a scaling plot of the mag- 
netization obtained with the multiple histogram technique according to its finite size scal- 
ing function (m w L»M (tL^) considering two possibilities: (i) the Ising-like behavior 
Rising _ 0.946(1), v = 1 and (3 = 0.125) and (ii) an order-disorder critical behavior with 
exponents v = 1.22(3) and /3 = 0.18(1) and critical temperature T c = 0.944(1). As can 
be seen, the scaling plot obtained assuming the Ising universality class does not describe 
our data as good as the results considering a new universality class. Besides, doing the 
same analysis with susceptibility and Binder's cumulant no significant deviations were ob- 
served between these two possibilities. Indeed, the values obtained in this study are in good 
agreement with those obtained for the same model in a bilayer system^ and for the dipolar 
Planar Rotator model^I. To clarify, in table [I] we show the exponents for the Ising model, 
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27 and 



the results obtained by Maier and Schwabl for the dPR model, the results of Refs. 
the results of this work. 

So far, everything corroborates to a order-disorder phase transition with non-conventional 
critical exponents. However, the scale relations^ a + 2/3 + 7 = 2 and vd = 2 — a are believe 
to be satisfied. Using the values shown in table [I] and the first relation we should have 
a = —0.51(7) and using the second relation a = —0.44(6) indicating the possibility that the 
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Figure 3: (Color online) Scaling plots of magnetization considering the Ising-like behavior (top) 
and an order-disorder transition characterized by the exponents shown in the last line of table [J 
(bottom). 



Model 



1 



Ising 2.269 1 1.75 0.125 (In) 

dPR (Maier) 1 1/2 -2 

AHd (bilayer) 0.890(4) 1.22(9) 2.1(2) 0.18(5) -0.55(15) 

dPR 1.201(1) 1.277(2) 2.218(5) 0.2065(4) -1.1(1) 

AHd (this work) 0.944(1) 1.22(3) 2.15(5) 0.18(1) -0.44(18) 

Table I: In this table we show the critical temperature and exponents for the 2D Ising model 3 ^ 
(first line), the results of Maier and Schwab!^ for the dPR model, the results of MC calculations 
in the bilayer AHd model with a cut-off in the interactions^, the results of MC calculations for 
the dPR model^ and the results of this work. 



specific heat does not diver ges. In deed, to have an better agreement between the results of 
this work and those of Refs. 
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J27jthe specific heat should be non-divergent. As one knows, 
to distinguish between an logarithmic divergence or an slowly power law divergence or even 
an non-divergent power law, many orders of magnitude are needed. Nevertheless, a careful 
analysis of the data could give us a clue. In figure [J] we show our data for the maxima of the 
specific heat as a function of the lattice size adjusted by two different methods. The dashed 
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Figure 4: (Color online) Specific heat maxima as a function of the lattice size. The dashed line is 
the best non-linear fit considering an logarithmic divergence and the dashed line shows the best 
non-linear fit considering a non-divergent power law behavior. 



line represents the best fit of a logarithmic divergence (aln(L) + b), the solid line is for a 
non-divergent power law behavior (—aL~ b + c). As can be clearly seen, the non-divergent 
power law describes better the data. Indeed, the x 2 /dof values obtained are 4.7 x 10~ 4 
for the logarithmic divergence and 1.4 x 10 -6 for the non-divergent power law. The value 
obtained for the exponent ajv from the adjust is —0.36(14), and it is also shown in tabled 
together with the results for the dPR and a bilayer AHd models. 



IV. DISCUSSION 



In this work we have studied the phase transition in the anisotropic Heisenberg model 
with dipolar interactions (AHd). Our main goal was to look for possible differences in 
the critical behavior of the planar to paramagnetic phase transition when full long-range 
dipolar interactions are considered or not. We have found that the use of the full long-range 



interaction by means of the Ewald summational^ lead to different results, i.e., whi 



introduction of a cut-off radius at five lattice spacings lead to a BKT transition (Ref. l2lh. the 



e the 



use of the Ewald summation lead to a order disorder transition with unusual exponents and 
a non divergent specific heat. Indeed, it would be interesting to present a detailed study of 
the effects in the critical behavior of the system as the cut-off radius increases. Nevertheless, 
this study is beyond the scope of this paper and will be addressed in a near future. On the 
other hand some points deserves a more detailed discussion. 



In what follows we present a discussion very similar to that done in Refs. 
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27. In 



the low temperature phase (T < T c ) the magnetization of the model does not display 
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any significant decrease as the lattice size is augmented. That is a clear indication of an 
order-disorder transition. This fact is in disagreement with the results obtained in the 
work of Rapini et al.— where was reported a BKT transition. Our results are very well 
described by a finite size scaling theory based on the existence of a low temperature phase 
with long range order and finite correlation length^. In a BKT phase transition there 
is no long-range order in the low temperature phase, consistent with the Mermin Wagner 
theorem^. Indeed, the results of Maleev— predict the existence of long-range order at low 
temperatures in the AHd model. Our results are consistent with this scenario. As discussed 
earlier, recent renormalization group results by Maier and Schwabl^ predicted that this 
system may belong to a new universality class, characterized by the presence of long-range 
order at low temperatures and by an exponential behavior of thermodynamic quantities in 
the vicinity of the critical temperature. In the Maier and Schwabl scenario the correlation 
length diverges as x ~ exp(6t 1 / 2 ) when the critical temperature (T c ) is approached from 
the high temperature side, similar to the behavior of the BKT phase transition, while the 
behavior of other thermodynamic quantities are given by powers of the correlation length. 
Nevertheless, our results for the AHd model are very well described by power law divergences 
of the thermodynamic quantities. As can be seen in figure [3], we have obtained a very good 
collapse of the curves from different lattice sizes for the magnetization and similar results 
were also found for the susceptibility and Binder's cumulant. These curves show that the 
critical exponents obtained and the conventional finite size scaling theory, that assumes a 
power law behavior of thermodynamic quantities, describe the Monte Carlo data accurately, 
indicating that the phase transition in the AHd model is a conventional order-disorder 
phenomenon with unusual critical exponents. In order to definitely rule out the possibility of 
this phase transition being in the new universality class proposed by Maier and Schwabl, we 
should make a comparison of our Monte Carlo results, using a finite size scaling theory based 
in their predictions and the conventional finite size scaling theory used here. Unfortunately, 
it is not very clear in the literature how to obtain a finite size scaling theory for exponential 
divergences. Using a simple replacement of the correlation length by the lattice size (in a 
manner similar to that made by Challa and Landau in reference^), which should be the first 
choice, does not give a good collapse of the curves, mainly because the determination of the 
critical temperature is quite imprecise in this case and the collapse of the curves depends 
appreciably on the value used for the critical temperature. In any case, using values for the 
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critical temperature close to the maxima of the susceptibility we were not able to obtain even 
a reasonable collapse of the curves. At this point one can argue: Why do renormalization 
group results not agree with Monte Carlo simulations? Actually, the RG study of Maier and 
Schwab!^ is based upon some approximations, for instance by using a continuous version of 
the model. Since the dipolar interactions have an intrinsic anisotropy, which depends in a 
complicated manner on the location of each spin in the lattice, the lattice geometry could 
have a strong effect in the system. The identification and a detailed discussion of the points 
of the RG study of the model that are behind the discrepancy between our results is beyond 
the scope of this paper. 

Although it was shown that the unusual exponents describes better the data, specially 
for the magnetization, the transition may be in the Ising universality class as well, since 
corrections to scaling were not taken into account and the lattice sizes used may not be 
large enough. Indeed, as can be seen in figure [3] the use of the Ising universality class 
exponents describes well the data for the largest lattices studied. Nevertheless, it does not 
seems to be a good choice simply disregard the data for the lattices with L = 20 and 40, 
leaving only two lattice sizes to be analyzed. Thus it is more prudent to not completely rule 
out the possibility of this phase transition to belong to the Ising universality class. On the 
other hand, our previous results for the same model in a bilayer system^ and the results 
for the dPR model^ were also well described by the same critical behavior found here, such 
that we still believe that this phase transition is more likely to belongs to a new universality 
class with unusual exponents. Studies in much larger lattices could remove this ambiguity, 
nevertheless the computational time needed for such a study turns it impracticable at the 
moment. 

Concerning the origin of the order disorder transition the question is even more compli- 
cated. The long-range order observed at low temperatures is expected to occur only when 
full long-range interactions are present. Nevertheless, in a recent study of the AHd in a 
bilayer system^ using a cut-off in the dipolar interaction, we found the same critical be- 
havior. The same results were also obtained in a study of the dipole planar rotator model 
in two dimensions^ with the dipole long-range interaction treated by means of the Ewald 
summation. In both cases the critical exponents agree quite well with those obtained in 
the present work (see tabled]). This observation indicates that the anisotropic character of 
dipolar interactions may be the main factor responsible for the observed critical phenom- 
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ena. Indeed, this observation is not new in the literature. As an example, Fernandez and 
Alonso^i stated that "Anisotropy has a deeper effect on the ordering of systems of classical 
dipoles in 2D than the range of dipolar interactions". In this work the authors found that 
the inclusion of a quadrupolar anisotropy drastically changes the phase transition behavior 
of a system of classical dipoles. Apparently, in our system the intrinsic anisotropy of dipo- 
lar interactions plays an essential role in the determination of the universality class of the 
AHd model. The possible new universality class is not surprising. In the theory of critical 
phenomenal it is expected that the critical exponents, and thus the universality classes, 
depend only on the spatial dimensionality of the system, the symmetry and dimensionality 
of the order parameter, and the range of the interactions within the system, characteristics 
not shared by the AHd model and models of well known universality classes. 

As a final remark we would like to stress that these results are much important when the 
critical behavior of magnetic systems with dipolar interactions is being considered. They 
show that the use of a cut-off radius in the evaluation of dipolar interactions may lead 
to erroneous results. This study may be a guide for future works in what concerns the 
introduction or not of a cut-off radius in the study of critical behavior of magnetic systems 
with dipolar interactions (see for example Ref. l43l ). 

We would like to thank Professor D P Landau for helpful discussions. Numerical calcula- 
tion was done on the Linux cluster at Laboratorio de Simulagao at Departamento de Ffsica 
UFMG. We are grateful to CNPq and Fapemig (Brazilian agencies) for financial support. 



* Electronic address: lucasmol@ufv.br 
' Electronic address: bvc@fisica.ufmg.br 



1 A. Berger, S. Mangin, J. McCord, O. Hellwig, and E. E. Fullerton, Phys. Rev. B 82, 104423 
(2010). 

2 S. D. Bader, Rev. Mod. Phys. 78, 1 (2006). 

3 D. P. Pappas, K.-P. Kamper, and H. Hopster, Phys. Rev. Lett. 64, 3179 (1990). 

4 R. Allenspach and A. Bischof, Phys. Rev. Lett. 69, 3385 (1992). 

5 L. A. M61, R. L. Silva, R. C. Silva, A. R. Pereira, W. A. Moura-Melo, and B. V. Costa, J. Appl. 
Phys. 106, 063913 (2009). 



12 



6 L. A. S. M61, W. A. Moura-Melo, and A. R. Pereira, Phys. Rev. B 82, 054434 (2010). 

7 E. Birsan, Journal of Superconductivity and Novel Magnetism 22, 711 (2009). 

8 S. Henning, F. Kdrmann, J. Kienert, W. Nolting, and S. Schwieger, Phys. Rev. B 75, 214401 
(2007). 

9 H. Leon and E. Estevez-Rams, Journal of Physics D: Applied Physics 42, 185012 (2009). 

10 E. Estevez-Rams and H. Leon, Journal of Magnetism and Magnetic Materials 322, 2874 (2010). 

11 D. Toscano, S. A. Leonel, R. A. Dias, P. Z. Coura, and B. V. Costa, J. Appl. Phys. 109, 076104 
(2011). 

12 D. Toscano, S. A. Leonel, R. A. Dias, P. Z. Coura, J. C. S. Rocha, and B. V. Costa, J. Appl. 
Phys. 109, 014301 (2011). 

13 V. L. Carvalho-Santos, W. A. Moura-Melo, and A. R. Pereira, J. Appl. Phys. 108, 094310 
(2010). 

14 F. A. Apolonio, W. A. Moura-Melo, F. P. Crisafuli, A. R. Pereira, and R. L. Silva, J. Appl. 
Phys. 106, 084320 (2009). 

15 K. De'Bell, A. B. Maclsaac, and J. P. Whitehead, Rev. Mod. Phys. 72, 225 (2000). 

16 C. Santamaria and H. T. Diep, J. Magn. Magn. Mater. 212, 23 (2000). 

17 M. Carubelli, O. V. Billoni, S. A. Pighm, S. A. Cannas, D. A. Stariolo, and F. A. Tamarit, 
Physical Review B (Condensed Matter and Materials Physics) 77, 134417 (pages 9) (2008). 

18 A. B. Maclsaac, K. De'Bell, and J. P. Whitehead, Phys. Rev. Lett. 80, 616 (1998). 

19 J. P. Whitehead, A. B. Maclsaac, and K. De'Bell, Physical Review B (Condensed Matter and 
Materials Physics) 77, 174415 (pages 18) (2008). 

20 H. Kwon, S. Hong, J. Seok, K. Bu, Y. Wu, Z. Qiu, and C. Won, Journal of Magnetism and 
Magnetic Materials 322, 2742 (2010). 

21 M. Rapini, R. A. Dias, and B. V. Costa, Phys. Rev. B 75, 014425 (2007). 

22 S. V. Maleev, Sov. Phys. JETP 43, 1240 (1976). 

23 P. Bruno, Phys. Rev. B 43, 6015 (1991). 

24 M. B. Taylor and B. L. Gyorffy, J. Phys.: Condens. Matter 5, 4527 (1993). 

25 P. G. Maier and F. Schwabl, Phys. Rev. B 70, 134430 (2004). 

26 L. A. S. M61 and B. V. Costa, Phys. Rev. B 79, 054404 (2009). 

27 L. A. S. M61 and B. V. Costa, J. Phys. Condens. Matter 22, 046005 (2010). 

28 S. K. Baek, P. Minnhagen, and B. J. Kim, Phys. Rev. B 83, 184409 (2011). 

13 



A. Carbognani, E. Rastelli, S. Regina, and A. Tassi, Phys. Rev. B 62, 1015 (2000). 

E. Rastelli, A. Carbognani, S. Regina, and A. Tassi, Eur. Phys. J. B 9, 641 (1999). 

J. F. Fernandez and J. J. Alonso, Phys. Rev. B 76, 014403 (2007). 

J.-J. Weis, Journal of Physics: Condensed Matter 15, S1471 (2003). 

Z. Wang and C. Holm, The Journal of Chemical Physics 115, 6351 (2001). 

R. P. Cowburn, D. K. Koltsov, A. O. Adeyeye, M. E. Welland, and D. M. Tricker, Phys. Rev. 

Lett. 83, 1042 (1999). 

J. C. S. Rocha, P. Z. Coura, S. A. Leonel, R. A. Dias, and B. V. Costa, J. Appl. Phys. 107, 
053903 (2010). 

R. F. Wang, C. Nisoli, R. S. Freitas, J. Li, W. McConville, B. J. Cooley, M. S. Lund, N. Samarth, 

C. Leighton, V. H. Crespi, et al., Nature 439, 303 (2006). 

D. P. Landau and K. Binder, A Guide to Monte Carlo Simulations in Statistical Physics (Cam- 
bridge University Press, New York, NY, USA, 2005), ISBN 0521842387. 

L. Onsager, Phys. Rev. 65, 117 (1944). 

V. Privman, ed., Finite Size Scaling and Numerical Simulation of Statistical Systems (World 
Scientific, 1990). 

N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 
M. S. S. Challa and D. P. Landau, Phys. Rev. B 33, 437 (1986). 

H. E. Stanley, Introduction to Phase Transition and Critical Phenomena (Clarendon Press - 
Oxford, 1971). 

L. A. S. M61, A. R. Pereira, and W. A. Moura-Melo, Physics Letters A 375, 2680 (2011). 



11 



